Computing Multi-Homogeneous Bezout Numbers is 



The mult i- homogeneous Bezout number is a bound for the number of solutions 
of a system of mult i- homogeneous polynomial equations, in a suitable product of 
projective spaces. 

Given an arbitrary, not necessarily multi-homogeneous system, one can ask for 
the optimal multi-homogenization that would minimize the Bezout number. 

In this paper, it is proved that the problem of computing, or even estimating 
the optimal multi-homogeneous Bezout number is actually NP-hard. 

In terms of approximation theory for combinatorial optimization, the problem 
of computing the best multi-homogeneous structure does not belong to APX, 
unless P = NP. 

Moreover, polynomial time algorithms for estimating the minimal multi-homo- 
geneous Bezout number up to a fixed factor cannot exist even in a randomized 
setting, unless BPP D NP. 

1 Introduction 

The multi-homogeneous Bezout number is a bound for the number of solutions of a 
system of multi- homogeneous polynomial equations. 

Estimating the number of isolated solutions of a polynomial system is useful for the 
design and analysis of homotopy algorithms |12| . Applications include problems in en- 
gineering like the design of certain mechanisms |15II18| or others, such as computational 
geometry. 

An application of multi-homogeneous Bezout bounds outside the realm of algebraic 
equation solving is discussed in 0], where the number of roots is used to bound geo- 
metrical quantities such as volume and curvature. 
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Abstract 



There is an important connection between root-counting and NP-completeness the- 
ory. Indeed, it is easy to reduce an NP-complete or NP-hard problem such as SAT, the 
Traveling Salesman problem, Integer Programming (and thus all other NP problems 
as well) to the question whether certain polynomial systems have a common zero. 

The best-known example giving an estimate for the number of roots of a polynomial 
equation is the Fundamental Theorem of Algebra. It was generalized to multivariate 
polynomial systems at the end of the 18th century by Etienne Bezout. The Bezout 
number bounds the number of (isolated) complex solutions of a polynomial / : C n i— * C n 
from above by the product of the degrees of the involved polynomials. However, in 
many cases this estimate is far from optimal. A well known example is given by the 
eigenvalue problem: Given a n x n matrix M, find the eigenpairs (A, it) eCxC™ such 
that Mu — Xu — 0. If we equate u n to 1, the classical Bezout number becomes 2 n ~ 1 , 
though of course only n solutions exist. 

The multi-homogeneous Bezout number provides a sharper bound on the number of 
isolated solutions of a system of equations, in a suitable product of projective spaces. 
The multi-homogeneous Bezout bound depends on the choice of a multi-homogeneous 
structure, that is of a partition of the variables (A, u) into several groups. 

In the eigenvalue example, the eigenvector u is defined up to a multiplicative con- 
stant, so it makes sense to define it as an element of P™ -1 . With respect to the eigenvec- 
tor A, we need to introduce a homogenizing variable. We therefore rewrite the equation 
as: XqMu — X\u = 0, and A = Xi/Xq. Now the pair (A : Ai) is an element of P 1 . The 
multi-homogeneous Bezout number for this system is precisely n. 

Better bounds on the root number are known, such as Kushnirenko's |l()j or Bern- 
stein's jSj- However, interest in computing the multi-homogeneous Bezout number 
stems from the fact that hardness results are known for those sharper bounds (see 
section l2~2l for details). 

Another reason of interest is that in many cases, a natural multi- homogeneous struc- 
ture is known or may be found with some additional human work. 

In this paper, we consider the following problem. Let n £ N and a finite icN" be 
given as input. Find the minimal multi-homogeneous Bezout number, among all choices 
of a multi-homogeneous structure for a polynomial system with support A: 



fi( z ) - J2 a eAfi<* z i 



. fn(z) — X^aeA fnatZ 1 1 Z 2 2 



(1) 



where the fi a are non-zero complex coefficients. 

Geometrically, this minimal Bezout number is an upper bound for the number of 
isolated roots of the system Q in C™. 

The main result in this paper (restated formally in section \2 . II below) is: 

Theorem 1. There cannot possibly exist a polynomial time algorithm to approximate 
the minimal multi-homogeneous Bezout number for up to any fixed factor, unless 
P = NP. 
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This means that computing or even approximating the minimal Bczout number up 
to a fixed factor is NP-hard. In terms of the hierarchy of approximation classes (see 
|2] and section ■ the minimal multi- homogeneous Bezout number does not belong 
to the class APX unless P = NP. 

Motivated by what is known on volume approximation (see section one could 
ask whether allowing for randomized algorithms would be of any improvement. 

Theorem 2. There cannot possibly exist a randomized polynomial time algorithm to 
approximate the minimal multi-homogeneous Bezout number for yj) up to any fixed 
factor, with probability of failure e < 1/4, unless BPP D NP . 

While the conjecture BPP ^ NP is less widely known outside the computer science 
community than the conjecture P ^ NP, its failure would imply the existence of 
probabilistic polynomial time algorithms for solving problems such as the factorization 
of large integers or the discrete logarithm. Most widespread cryptographic schemes are 
based on the assumption that those two problems are hard. 



2 Background and Statement of Main Results 

2.1 Bezout numbers 

In the definition of Q, we assumed for simplicity that each equation had the same sup- 
port A. In general, a system f{z) of n polynomial equations with support (A\, . . . , A n ) 
is a system of the form: 



H Z ) = HaeA.hocZ^Z 



"l an 
' ' ' *« 

(2) 



fn(z) — YlaeA n fna.Z l 1 Z 2 2 • • • z" n , 

where the coefficients fi a are non-zero complex numbers. 

A multi-homogeneous structure is given by a partition of {1, . . . , n} into (say) k sets 
I\, . . . , Ik. Then for each set Ij, we consider the group of variables Zj = {zi : i € Ij}. 

The degree of fi in the group of variables Zj is 



, dcf \ ^ 

da = max > ct] 



When for some j, for all i, the maximum is attained for all a e ^, we say 
that J5J) is homogeneous in the variables Zj. The dimension of the projective space 
associated to Zj is: 

dcf J #Jj — 1 if (0 is homogeneous in Zj, and 
3 I ^I'j otherwise. 

We assume that n = X^=i a i- Otherwise, we would have an undetermined (n < 

Y]j—i <ij) or overdetermined (n > Y)j—i &j) polynomial system, and multi- homogeneous 
Bczout numbers would have no meaning. 
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The multi-homogeneous Bezout number Bez(Ai, . . . , A n ; I\, . . . , Ik) is the coefficient 
of Ylj—i Cj J m the formal expression 

11!' i Y., i ( see lEllISlini). It bounds the 
maximal number of isolated roots of in P ai x ■ ■ • x P afc . Therefore it also bounds 
the number of finite roots of @, i.e. the roots in C™. 

In the particular case where A = A\ = ■ ■ ■ = A n there is a simpler expression for the 

multi-homogeneous Bezout number Bez(A; I±, . . . , Ik) == f Bez(Ai, . . . , A n ; I\, . . . , Ik), 
namely: 

where ci, = c?y (equal for each i) and the multinomial coefficient 

/ n \ def ?i! 

\ «i "2 • ■ ■ fli / a\\ a%\ ■■■ cik ! 

is the coefficient of n* = i C" fc in (Ci + • ■ • + (k) n (recall that n = a^). 

Heuristics for computing a suitable multi-homogeneous structure (ii,... , Ik) given 
Ai, . . . ,A n are discussed in [TTlfTT^j . Surprisingly enough, there seems to be no theoret- 
ical results available on the complexity of computing the minimal Bezout number. It 
was conjectured in p. 78] that computing the minimal multi-homogeneous Bezout 
number is NP-hard. 

Even, no polynomial time algorithm for computing the multi- homogeneous Bezout 
number given a multi-homogeneous structure seems to be known (see |13l p.240]). 

This is why in this paper, we restrict ourselves to the case A = A\ = ■ ■ ■ = A n . This 
is a particular subset of the general case, and any hardness result for this particular 
subset implies the same hardness result in the general case. 

More formally, we adopt the Turing model of computation and we consider the 
function: 

Bez :n,k,A,h,...,I k ^> Bez(A;/i, ...,Ik) 

where all integer numbers are in binary representation, and A is a list of 7i-tuples 
(ai, . . . , a n ), and each Ij is a list of its elements. In particular, the input size is bounded 
below by n^Ai and by max Qj i |~log 2 ai\ . Therefore, Bez(A; I\ , • • • , Ik) can be computed 
in polynomial time by a straight-forward application of formula 0. As a matter of 
fact, it can be computed in time polynomial in the size of A. 

Problem 1 (Discrete optimization problem). Given n and A, compute 

minBez(A; I) , 

where I = (7i, . . . , Ik) ranges over all the partitions of {1, ... ,n}. 

Problem 2 (Approximation problem). Let C > 1 be fixed. Given n and A, 
compute some B such that 

BC^ 1 < minBez(A;/) < BC 
Again, I = (ii, . . . , Ik) ranges over all the partitions of {1, . . . , n}. 
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In the problems above, we are not asking for the actual partition. 

Theorem 1 (restated). Problem^ is NP-hard. 

This is actually stronger than the conjecture by Li and Bai JJ, that corresponds 
to the following immediate corollary: 

Corollary 1. Problem^ is NP-hard. 

2.2 Other bounds for the number of roots 

Kushnirenko's Theorem ^U] bounds the number of isolated solutions of in (C*) n by 
n\ Vol Convj4, where Convj4 is the smallest convex polytope containing all the points 
of A. 

This bound is sharper than the Bezout bound, but the known hardness results are 
far more dramatic: In Khachiyan proved that computing the volume of a polytope 
given by a set of vertices is #P-hard. 

There is a large literature on algorithms for approximating the volume of a con- 
vex body given by a separation oracle. The problem of approximating the volume of 
a polytope in vertex representation can be reduced to the latter by standard linear 
programming techniques. 

It is known that no deterministic algorithm can approximate the volume in polyno- 
mial time f|14j). However, randomized polynomial time algorithms are known for the 
same problem [71H^]. 

The same situation seems to be the case regarding the estimation of the mixed 
volume jHj, which gives the actual number of solutions in (C*) n for generic polynomials 
of the form © 0. 

2.3 Probabilistic algorithms 

A probabilistic machine is a machine that has access to random bits of information, each 
random bit costing one unit of time. Each random bit is an independent, uniformly 
distributed random variable in {0, 1}. In that sense, a probabilistic machine is a machine 
that flips a fair coin, as many times as necessary, spending one unit of time at each flip. 

We can therefore speak of the probability that the machine returns a correct result. 

The class BPP is the class of decision problems (X, X ycs ) such that there is a 
probabilistic machine and a constant e < 1/2 that will: 

(i) Decide in polynomial time if x € X. 

(ii) Output YES or NO, in polynomial time. 

(iii) For every x, the output is the correct answer to the question: does x S X yes ? 
with probability > 1 — e. 

Notice that we can improve the probability that the result is correct by running the 
same machine several times. Therefore, in the definition above, we may as well take 
e = 1/4. 

More generally, a probabilistic machine solves a certain problem (e.g. Problem^ in 
polynomial time with probability > 1 — e if and only if it always terminates in polynomial 
time, and the answer is correct with probability 1 — e. 
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Theorem 2 (restated). There is no e < 1/2 and no probabilistic machine solving 
Problem\^ with probability 1 — e, unless BPP D NP. 

2.4 Approximation classes 

A theory of complexity classes appropriate for the study of combinatorial optimization 
problems is described in [2] . Problem ^ fits naturally in the class of combinatorial 
optimization problems. In this context, Problem ^ is characterized by: 

(i) A set of instances, given by the set of pairs (n, A), n G N, Ac N n finite and 
non-empty. 

(ii) For every instance (n, A), a set of feasible solutions, namely the set of partitions 
I=(h,...,I k ) of{l,. ..,n}. 

(hi) An objective function (to minimize), Bez(A; I). 

The class NPO of combinatorial optimization problems is analogous to the class 
NP of decision problems. Problem ^ belongs to that class: 

(1) The size of each feasible solution is polynomially bounded on the size of each 
instance. 

(2) Given an instance (n, A) and a string w, it can be decided in time polynomial in 
(n, A) whether w encodes a feasible solution I = (Ii, . . . , Ifc). 

(3) The objective function can be computed in polynomial time. 

The class APX of approximable problems in NPO is defined as the subset of NPO 
for which there is some C > 1 and a polynomial time algorithm such that, given an 
instance of the problem (say n, A)) produces a feasible solution I such that the objective 
function applied to that solution approximates the minimum up to a factor of C . 

Theorem [5] admits as a corollary: 

Corollary 2. Problem^ does not belong to APX, unless P = NP. 

Our result actually holds even if we do not require the algorithms to compute a 
feasible solution. 

3 Proof of the Main Theorems 

3.1 Prom graph theory to systems of equations. 

Definition 1. A k-coloring of a graph G = (V, E) is a partition of the set of vertices 
V into k disjoint subsets ("colors") Ij, so that adjacent vertices do not belong to the 
same "color" Ij. 

Problem 3 (Graph 3-Coloring). Given a graph G = (V,E), decide if there exists a 
3-coloring of G. 



G 




Figure 1: In this example, A{H) = {(0, 0, 0, 0), (1, 0, 0, 0), (0, 1, 0, 0), (0,0,1,0), 
(0,0,0,1), (1,1,0,0), (1,0,1,0), (0,1,1,0), (0,0,1,1), (1,1,1,0)}. A possible polyno- 
mial with that support would be 1+V\ + V2 + «3 + V4 + V1V2 + V1V3 + V2V3 + V3V4 + V1V2V3. 



It is known since [B] that the Graph 3-Coloring Problem is NP-hard (see also 
We will actually need to consider an equivalent formulation of the Graph 3-coloring 
problem. 

Recall that the cartesian product of two graphs G\ = (Vi, Ei) and G2 — (V'2, E2) is 
the graph G\ x G2 = (V"i x V~2,E) with V2), (v[, v' 2 )) & E if and only if v\ = v[ and 
(V2,v' 2 ) e E 2 or v 2 — v' 2 and (vi,v[) € E\. 

Also, let K 3 denote the complete graph with 3 vertices. 

Lemma 1. The graph G admits a 3-coloring if and only if the graph G x if 3 admits a 
3-coloring I = {h,h,h) with #7j = #/ 2 = #J 3 = \G\. ' 

Proof. G admits a 3-coloring if and only if G x K3 admits a 3-coloring. Moreover, any 
coloring J of G x if 3 satisfies #ii = #1% = #13. □ 

To each graph H = (V,E) we will associate two spaces of polynomial systems. 
Each of those spaces is characterized by a support set A = A(H) (resp. A(H) 1 ) to be 
constructed and corresponds to the space of polynomials of the form with complex 
coefficients. Of particular interest will be graphs of the form H = G x K3. 

We start by identifying the set V of vertices of H to the set {1, . . . ,m}. Let K s 
denote the complete graph of size s, i.e. the graph with s vertices all of them pairwise 
connected by edges. 

To each copy of K s , s = 0, . . . , 3 that can be embedded as a subgraph of H (say the 
subgraph generated by {v\, • ■ ■ ,v s }) we associate the monomial 

(the empty graph K corresponds to the constant monomial). Then we consider the 
linear space generated by all those monomials (Figure^). Therefore, the support A(H) 
is the set of all e Vl + ■ ■ ■ + e Vs c W n such that < s < 3 and . . . , v s } induces a 
copy of K s as a subgraph of H . Here, e, denotes the i-th vector of the canonical basis 
of W 1 , 

Given a set A, we denote by A 1 the /-fold i cartesian product of A. 
The two spaces of polynomial systems associated to a graph H will be the polynomial 
systems with support A(H) and A{H) 1 . 
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Remark that none of the two classes of systems above is homogeneous in any possible 
group of variables (because we introduced a constant monomial). Therefore, in the 
calculation of the Bezout number for a partition J, we can set aj = jfclj. 

Lemma 2. Let I be fixed. Then, there is a polynomial time algorithm to compute A(H) 
and A(H) 1 , given H . 

3.2 A gap between Bezout numbers 

In case the graph H admits a 3-coloring J = (Ji, I 2 , I3), any corresponding polynomial 
system is always trilinear (linear in each set of variables) . If moreover H is of the form 
H = G x K3 with \G\ = n. the cardinality of the Ij is always n, and formula 
becomes: 

Bez(A(GxK 3 );I)= ( 3n ) 
v v ' ' V n n n J 

The crucial step in the proof of Theorem ^ is to show that 

Bez(A(GxK 3 );I)>~( 3n ) 
x v ' ' 3 V n n n J 

unless k = 3 and J is a 3-coloring of G x K3. 

In order to do that, we introduce the following cleaner abstraction for the Bezout 
number: if k € N and a = (aj., . . . , a^) € N k are such that Sj=i a j = we se ^ 

dcf/ 3n \ Jr a, 

\ ax a 2 ■ ■ ■ ak / 1L I n 

\ ' j=i 

Lemma 3. If H = G x K3 and I — (Ji, . . . ,/&) is a partition of the set {1, . . . ,3n} of 
vertices of H , then 

Bez(A(H);I) > B(a) 

with aj = #Ij . 

Proof. Consider the n disjoint copies of K3 in H = G x K3 induced by the nodes of G. 
By the pigeonhole principle, there is at least one of those copies with at least |~aj/n] 
elements of Ij . Hence, the degree dj in the j-th group of variables is at least \aj/n \ . □ 

The main step towards establishing the "gap" is the following Proposition: 

Proposition 1. Let n, k 6 N n and let ai > a 2 > ■ ■ ■ > ak > 1 be such that 53,- =1 dj — 
3n. Then, either k — 3 and a\ = a 2 = 03 = n, or: 

B(a) > ^B{n,n,n) . 

Moreover, this bound is sharp. 

The proof of Proposition 2] is postponed to section 0] 
Putting it all together, 
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Lemma 4. Let G be a graph and n 



\G\. If G admits a 3-coloring, then 



Otherwise, 




) 



Proof. According to Lemma ^ G admits a 3-coloring if and only if G x K3 admits a 
3-coloring. 

If I = (ii, 12,13) is a 3-coloring of G x K3, then 



If I = (h, . . . , Ik) is not a 3-coloring of G x K3, then we distinguish two cases. 
We set aj = #Ij . 

Case 1: a = (n,n,n) and hence k — 3. Then the degree in at least one group of 
variables is > 2, and 



where the first inequality follows from Lemma Eland the second from Proposition 
In both cases, 



Lemma 0] would be sufficient to prove a weaker version of Theorem ^ where the 
factor C in problem |3 is less than 4/3. 

3.3 Improving the gap 

In order to obtain a proof valid for any C the idea is to increase the gap by considering 
several copies of a polynomial system, but each copy in a new set of variables. This 
idea works out because of the special multiplicative structure of the multi-homogeneous 
Bezout number. We will need: 

Proposition 2. Let to, I <E N. Let A C N m be finite and assume that G A. Then, 





Case 2: a / (n,n,n). Then 





□ 
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Proof. 1. Let I = (7i, • • • , Ik) be the partition of {1, ... , to} where the minimal Bezout 
number for A is attained. 

This induces a partition J — (Jj S )i<j<k,i<s<i of {1, . . . , to} x {1, . . . , I}, given by 
Jj S = Ij x {s}. Identifying each pair (i, s) with i + tos, the Jj S are also a partition of 
{1, . . . ,lm}. 

By construction of A 1 , the degree dj S in the variables corresponding to Jj s is equal 
to the degree dj of the variables Ij in A. 

The systems corresponding to A and A 1 cannot be homogeneous for any partition, 
since <G A and € A 1 . Then we have dj = #Ij = a JS for any s. Therefore, 



mmB6z(A l ,K) < Bez(A l ,J) 




a\ a 2 ■ ■ ■ afc 



IK 



^minBez(A; 1)^ 



2. Now, suppose that the minimal Bezout number for A 1 is attained for a partition 
J = ( Ji, • • • , J r ). We claim that each J t fits into exactly one of the I sets {1, . . . , to} x 
{*}. 

Suppose this is not the case. Assume without loss of generality that Ji splits into 
K C {1, . . . ,to} x {1} and L C {1, . . . ,to} x {2, . . . ,^}, both if and L non-empty. 

If d K denotes the degree in the if-variables and d^ the degree in the L variables, 
then di = d K + d L . Also, a\ = a K + a L where a K is the size of K and is the 
size of L. The multi-homogeneous Bezout number corresponding to the partition J' = 
(K,L,J 2 ,--- ,Jr) is: 

Bcz(A'; J') = ( Zlm )d a Jd a S fj d a / 

V ' \ a K a L a 2 ■ ■ ■ a r J K L 11 i 



Therefore, 

B6z(A l ;j') \ a 



ai ^d a K K d a L L 



Bez(A l ,J) " {d K + d L ) ai 



< 1 



and the Bezout number was not minimal, thus establishing the claim. 

3. Denote by J = U l s=1 J^ the partition minimizing the Bezout number correspond- 
ing to A 1 . In the notation above, we assume that is a partition of {1, ... , to} x {s}. 
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In that case, 

Bez(A ! ; J) . („ m te, .. m )S((^"^)n^) 

= ( )riBez(AJW) 
\ m m ■ ■ ■ m 

v 7 s=l 

> ( lm VminBez(A;I)V 
\ m m ■ ■ ■ m J \ i J 

□ 

Combining Lemma 0] and Proposition we established that: 
Lemma 5. Lei G be a graph and n — \G\. Let i £ N. If G admits a 3-coloring, then 

mmBez(A(GxK 3 ) l ,J)=( 3nl )( 3n ) 
j \ on on ■ ■ ■ on J \ n n n I 

Otherwise, 

^B6,(A(G,K 3 y,J)> (!)'( 3n3 ™. 3n )( n 3 : n )' 

Proof of Theorem^ Assume that ApproxBez is a deterministic, polynomial time algo- 
rithm for solving problem [21 i.e., for estimating the Bezout number up to a factor of 
C. 

Then the following algorithm decides Graph 3-coloring (Problem |3J) in polynomial 
time: 

Algorithm 1 (Decides Graph 3-coloring problem). 
Input : a graph G of size n. 

Output : YES if G admits a 3-coloring, NO otherwise. 
Constants: 1= 2 }ggf /3 . 

1. Compute 

ApproxBez (A(G x K 3 ) 1 ) 

P '~ ( onl U 3n V 
\ 3n 3n ■ ■ ■ 3n J \ n n n J 

2. If p 2 < C then Output YES, else Output NO. 

By our choice of the constant I, \[C < (4/3)'. Therefore, Lemma |3] asserts that the 
output of algorithm n is correct. 

The bit-size of the numbers that occur when computing the denominator of line 
2 are bounded above by 0(3nl log(3nZ)). The size of the graph G x Kz is 0(n), and 
Lemma |2 says that A 1 can be computed in polynomial time. 

It follows that Algorithm ^ runs in polynomial time. Since Graph 3-coloring is 
NP-complete, we deduce that P = NP. □ 
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Proof of Theorem^ Assume now that ApproxBez is a probabilistic polynomial time 
algorithm for solving problem [21 which returns a correct result with probability 1 — e, 
e < 1/4. 

Then Algorithm ^ will return the correct answer for the Graph 3-coloring Problem, 
with probability at least 1 — e. This implies that Problem[3]is actually in BPP. □ 



4 Proof of Proposition [T] 

We will need the following trivial Lemma in the proof of Proposition ^ 
Lemma 6. Let x, n e N. Then. 

> 1 + ((n — x) mod n) . 

In 'particular, the left-hand side is > 2 whenever n }(x, and is always > 1. 
Proof. Since n [~~] = x + (n — x) mod n, we have: 

(n — x) mod n\ x 







n 


a; / 



> 1 + ((n — x) mod n) 



Also, we will make use of the Stirling Formula (6.1.38)]: 

/ -1-1 | 0(aQ 

where < 8(x) < 1. 

Proof of Proposition^ The ratio between -B(a) and i?(?T,, n, n) is: 



B(a) 



-B(ri, n, n) 



n 



)t «.! rjt 



ai! a 2 ! • • • a k \ 



From Stirling formula (0J it follows immediately that: 



S(a) 



Bin, n, n) 



'2n 



-3-k 



n 



-371+4 



2-i 12 



□ 



(4) 



(5) 



Now we distinguish the cases k = 1, fc = 2, and fc > 3. The first two cases are easy: 
Case 1: If fe = 1, then ai = 3n and JSJ) becomes: 

S(a) „ n sm f(3r») 
— = 2-7T — = e 4 " 

B{n,n,n) V3 

which is bounded below by ^e" 1 / 36 ~ 3.528218766. 
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Case 2: If k = 2, Lemma [S] implies that 



B(n,n,n) ^/ol 02 

Since ^/oTo^ < = TP we obtain: 

„, B(a) r > IV^e- 1 / 6 ~ 1.414543350 



Case 3: Let k > 3. If 03 = n, then k — 3 and ai = 02 = &3 = n>, so there is nothing 
to prove. Therefore, we assume from now on that 03 < n. 

We separate the right-hand side of (JSJ into two products, the first for j — 1,2,3 and 
the second for j > 4. Equation JSJ becomes now: 



B(a) 



B(n, n, n) 



nil? 



n 2 



8(n) _ v a e <°j> 
^ 4?i ^? = 1 12a ^ 



(6) 



2tt 



3-fc 



11 aj +i 



= 2^ 3= 4 12a,,- 



J=4 ab- 



using the fact that a, < n for j > 4. In case fc = 3, the second factor in equation JJjJ 
above is equal to one. 

Since 03 < n, n )( 0,3 and Lemma [5] implies that for 03 < n 



nil? 

3=1 



— I >2 



Moreover, ^/aia^ai < (a\ + 0,2 + 0,3) /3 < n, so the first factor of the right-hand side 
of (JIU can be bounded below by 



k 

n 

3=1 



— e 4 " ^ 3=1 12a o 



> 2e 



-1/4 



1.557601566 



If k = 3 we are done. Otherwise, we notice that since the are non-increasing, 
< for all j > 4. In order to bound the second factor of iJBJ, we will need the 



following technical Lemma: 

Lemma 7. Let n, x £ N ar, 

(6,4), (7,5), (8,6)}, we have: 



Lemma 7. Let n, x S N and let x< ^f. Then, unless (n,x) € {(2, 1), (3, 2), (4, 3), 



27T£ x+ t 



-e 12- > 1 



(Proof is postponed). 
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n 




3n 


a 


j 


5(a) 


B(n, n, n) 


aw 

B(n,n,n) 


2 


1 


6 


111111 
2 1111 

2 2 11 

3 111 


4,5,6 
4,5 
4 
4 


720 
360 
180 
120 


90 


8 
4 
2 

4 


3 


2 


9 


2 2 2 2 1 

3 2 2 2 


4 
4 


22680 
7560 


1680 


— k — 


4 


3 


12 


3 3 3 3 


4 


369600 


34650 


— h — 

3 


6 


4 


18 


4 4 4 4 11 

4 4 4 4 2 

5 4 4 4 1 

5 5 4 4 

6 4 4 4 


4 
4 
4 
4 
4 


19297278000 
9648639000 
3859455600 
771891120 
643242600 


17153136 


1125 

1125 
2 

225 
45 

75 


7 


5 


21 


5 5 5 5 1 

6 5 5 5 


4 
4 


246387645504 
41064607584 


399072960 


30%7 

5 

1029 


8 


6 


24 


6 6 6 6 


4 


2308743493056 


9465511770 


45 



Table 1: Ratios for all the exceptional pairs (n,a). 



Therefore, unless some of the pairs (n,a,j), j > 4 belong to the exceptional subset 
{(2, 1), (3, 2), (4, 3), (6, 4), (7, 5), (8, 6)}, we have: 

5(a) i 
v ' > 2e~2 ~ 1.557601566 . 



B(n, n, n) 



Finally, we consider the values of n and a where some (n,aj), j > 4, is in the 
exceptional subset. All the possible values of n and a are listed in tabled The ratio is 
always > 4/3, and the value of 4/3 is attained for n — 2 and a = (3, 1, 1, 1). □ 

Proof of Lemma^ Let 



9n{x) = log( ^^T e " 12a 



n 



— x log n — x log x log x log 2tt 

5 6 2 6 12x 2 5 

(see figure We first consider values of x > 7. By hypothesis, n/x > 4/3 so 
logn — log a; > log(4/3), and therefore g n {%) > h(x), where: 

h(x) = xlog(4/3) - i log a; - - \ log27r 
2 12a; I 

(see Figure |21 also) . Notice that h(x) is independent of n. The derivative of h is 
The numerator vanishes at 



l±lv/l-4/31og(4/3) 
41og(4/3) 
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Figure 2: Plots of g n {x) and h(x). 



Numerically, x ~ 0.1867281114 on~ 1.551301638. Therefore, the function h(x) is 
increasing for x > 2. Again, numerically h(7) ~ 0.1099761345 and therefore, if x > 7 
we always have: 

e ff»(aO > 1.H62 > l 
Now we consider x < 6. Having g n > 0, is equivalent to: 
n > n Q (x) = xe^ losx+ Tt?+^ l °s^ 

At this point, we proved that g n (x) is positive, except possibly for pairs (n, x) with 
1 < x < 6 and 3X < n < no(-' : ' ; )- The values of uq are tabulated in Table [5] From 
Table [3 it is clear that the only exceptions are those listed in the hypothesis. 



X 




n (x) 


Possible n's 


1 


1.333333333 


2.724464424 


2 


2 


2.666666666 


3.844857634 


3 


3 


4 


4.939610298 


4 


4 


5.333333333 


6.016610872 


5 


5 


6.666666666 


7.081620345 


6 


6 


8 


8.137996302 


8 



Table 2: Possible values of n for x small 



□ 
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